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Studies of first-order phase transitions through the use of the exact renormalization group 
CM are reviewed. In the first part the emphasis is on universal aspects: We discuss the uni- 

versal critical behaviour near weakly first-order phase transitions for a three-dimensional 
^ ■ model of two coupled scalar fields - the cubic anisotropy model. In the second part we 

Q review the application of the exact renormalization group to the calculation of bubble- 
nuclcation rates. More specifically, we concentrate on the pre-exponential factor. We 
, discuss the reliability of homogeneous nucleation theory that employs a saddle-point 

expansion around the critical bubble for the calculation of the nucleation rate. 

. 1 Weakly first-order phase transitions 

i> , 

Weakly first-order phase transitions appear in many statistical systems (such as su- 
perconductors or anisotropic systems) and in high-temperature field theories. Many 
cosmological phase transitions, such as the electroweakEl, fall in this category. How- 
ever, a reliable quantitative description of these phenomena is still lacking. 

A simple example of a system with an arbitrarily weakly first-order phase tran- 
sition is the cubic anisotropy modelBa. In field-theoretical language, it corresponds 
to a theory of two real scalar fields <p a (a = 1,2) in three dimensions, invariant 
O ' under the discrete symmetry (1 «-» — 1, 2 <-> —2, 1 <-> 2). It can also be considered as 

the effective description of the four-dimensional high-temperature theory with the 
same symmetry, at energy scales smaller than the temperatureu. 

The properties of the weakly first-order phase transitions in this modeLhave 
been studied within the e-expansion El or through Monte Carlo simulations In 
particular, universal amplitudes have been computed, which describe the relative 
discontinuity of various physical quantities along the phase transition, in the limit 
when the transition becomes arbitrarily weakly first order. A discrepancy has been 
observed in the predictions for the universal ratio of susceptibilities x+/X— on cither 
side of the. phase transition, obtained through Monte Carlo simulations or the e- 
expansionQ. The Monte Carlo simulations predict x+lx- = 4-1(5), while the first 
three orders of the e-expansion give X+/X- = 2.0, 2.9, 2.3 respectively. 

We summarize here the results of an alternative approach a that employs the 
exact renormalization group. It is based on the effective average action u, which 
is a coarse-grained free energy with an infrared cutoff. More precisely, Tk incor- 
porates the effects of all fluctuations with momenta q 2 > k 2 , but not those with 
q 2 < k 2 . In the limit k — > 0, the effective average action becomes the standard 
effective action (the generating functional of the 1PI Green functions), while at a 
high momentum scale (of the order of the ultraviolet cutoff) k — A — > oo, it equals 
the classical or bare action. An exact non-perturbative flow equation determines 
the scale dependence of IV 
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The flow equation is a functional differential equation, and an approximate 
solution requires a truncation Our truncation is the lowest order in a systematic 
derivative expansion of FfcEiEJ 

r fe = J ePxlu h (j>i,p2) + \z k (Wfadntf + 0"M^ 2 )} . (1) 

Here p a = ^4> a 4> a & n( A the potential Uk(pi, P2) is symmetric under the interchange 
1 <-> 2. The wave- function renormalization is approximated by one fc-dependent 
parameter The truncation of the higher derivative terms in the action is ex- 
pected to generate an uncertainty of the order of the anomalous dimension 77. For 
the model we are considering, 77 ~ 0.035 and the induced error is small. 

The fixed-point structure of the theory is more easily identified if we use the 
dimensionless renormalized parameters 

p a = Z k k~ 1 p a u k (pi,p 2 )= k~ 3 U k (pi,p 2 )- (2) 

The evolution equation for the potential can now be written in the scale- independent 
forma 

d 11 

■j£Uk(pi,fc) = -3u fe + (l + 7j)(/5iui + p 2 u 2 ) - ^Lg(m^) - ^^(n^), (3) 

where t = ln(fe/A). The anomalous dimension r\ is defined as dh\Z k /dt = —77. 
For the part o£,the phase diagram of interest, it is constant, 77 ~ 0.035, to a good 
approximationtl The quantities rh\ 2 are the eigenvalues of the rescaled mass matrix 
at the point (p±, P2) 



2fh\ 2 = ui + U2 + 2uiipi + 2U22P2 ± y (ui - u 2 + 2u\\p\ - 2U22P2) 2 + ld>u\ 2 pip2. 

We use the notation u\ = duk/dpi, u\2 — & 1 u k jdp\dp2 1 etc. The function L\{w) is 
negative for all values of w. Also is monotonically decreasing for increasing 

w and introduces a threshold behaviour in the evolution. For large values of rh 2 the 
last two terms in eq. (^) vanish and the evolution of Uk stops Ui3. 

The initial condition for the integration is provided by the bare potential, which 
is identified with the effective average potential at a very high scale k = A. We use 
a bare potential of the form 

"a(/5i,/3 2 ) = ^A A - ka) 2 + (P2 - «a) 2 | + (1 + x A )\ A pip~2 (5) 

The phase diagram of the theory §i@ has three fixed points that govern the 
dynamics of second-order phase transitions. They are located on the critical surface 
separating the phase with symmetry breaking from the symmetric one. The most 
stable of them corresponds to a system with an increased 0(2) symmetry. It can be 
approached directly from a bare action given by eq. (Q) with x\ — 0. The other two 
are Wilson-Fisher fixed points, corresponding to two disconnected .^-symmetric 
theories. One of them can be approached from a bare action with xa = — 1, while 
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Figure 1: The derivative of the rescaled potential along the pi axis, as the coarse-graining scale 
is lowered. 

the second requires x\ = 2 ^|. Flows that start with — 1 < < 2 and near the 
critical surface eventually lead to the 0(2)-symmetric fixed point. For x\ > 2 or 
x\ < — 1 the evolution leads to a region of first-order phase transitions. If x\ is 
chosen slightly larger than 2 or slightly smaller than —1 the phase transitions are 
weakly first order. The evolution first approaches one of the fixed points before a 
second minimum appears in the potential. 

te partial differential equation (^) with r\ — 0.035 can be integrated numeri- 
We are interested in the region x.> 2. The region x < — 1 can be mapped 
on it through a redefinition of the fields □. We concentrate on the potential along 
the p\ axis and we define Vk{pi) = Uk(pi, f>2 = 0). In fig. 1 we present the evolution 
of vi — dvk/dpi, starting at a very high scale k = A with a bare potential given by 
eq. (||). All dimensionful quantities are normalized with respect to A. The initial 
coupling Aa is chosen arbitrarily, while the minimum k\ of the bare potential is 
taken very close to the critical value n cr that separates the phase with symmetry 
breaking from the symmetric one. For \6n\\ = \ka — « cr | -C 1 the system spends 
a long "time" t of its evolution on the critical surface separating the two phases. 
The initial value x\ is taken slightly larger than the fixed-point value 2. After the 
initial evolution (dotted lines) the potential settles down near the Wilson-Fisher 
fixed point (solid lines). This fixed point, that has a mixing term between the two 
fields with x — 2, is repulsive in the x direction. Eventually x evolves towards 
larger values. This forces the potential to move away from its scale-independent 



a A redefinition of Ahe fields demonstrates that this choice corresponds to two disconnected Z2- 
symmetric theories □. 
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Figure 2: The potential along the pi axis, as the coarse-graining scale is lowered. 



form (dashed lines). At some point in the subsequent evolution the curvature of the 
potential at the origin vi(pi = 0) becomes positive. This signals the appearance 
of a new minimum there, and the presence of a radiatively-induced (or fluctuation- 
driven) first-order phase transition. The evolution of the potential after it moves 
away from the fixed point is insensitive to the details of the bare potential. It is 
uniquely determined hx \Sk\\ = \ka~ n cr \ <§; 1 and Sxa — x\ — 2 <C 1, and displays 
universal behaviour 113. 

In the limit k — > and in the convex regions, the rescaled potential Uk grows 
and eventually diverges in such a way that Uk becomes asymptotically constant, 
equal to the effective potential U = Uq. This is apparent in fig. 2, where the 
potential along the p\ axis is plotted. The evolution of the non-convex part of the 
potential (between the two minima) is related to the issue of_the convexity of the 
effective potential. This part should become flat for k — ► 0c3. The approach to 
convexity is apparent in fig. 2, even though we have not followed the evolution all 
the way to k = 0. 

The relative magnitude of \6k\\ — \ka — n cr \ -C 1 and Sx a = x\ — 2 <C 1 results 
in different types of evolution. For the type of behaviour depicted in figs. 1 and 2 
one must take \5k\\ <C 5x\. In the opposite limit, \8k\\ 3> 5x\, the system leaves 
the critical surface before the mixing x between the fields evolves away from the 
fixed-point value 2. As a result, a second minimum never appears at the origin. 
Instead, the only minimum of the effective potential Uq is located either at zero 
(symmetric phase) or away from it (phase with symmetry breaking), depending 
on the sign of Sn\. The resulting phase transition is second order and occurs for 
Ska = 0. 
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We are interested in the universal ratio of susceptibilities x+/X- on either 
side of the phase transition. This ratio depends on the value of \Ska\/Sxa- For 
\Ska\/Sxa > 1 we have x ~ 2 during the whole evolution. The universal quantities 
characterizing the second-order phase transition are determined by the Wilson- 
Fisher fixed point. We calculate x+/X- by integrating the evolution equation and 
evaluating cPUo/d<jy( — x 1 at the minimum, for Ska — T e with f«l. We obtain 
X+/x- = 4.1. The evaluation of the same quantity through the e-expansion or 
an expansion at fixed dimension gives X+/X- = 4.8(3) Ej, whereas experimental 
information gives \+/X- = 4.3(3) E3. The difference with the results of our method 
is due to the omission of higher derivative terms in the effective average action of 
eq. (fil) and approximations that have been used for the numerical integration of 
eq. (|p. 

For \8ka\I 8xa <C 1 the potential develops a second minimum at the origin 
during the later stages of the evolution. The phase transition is approached by 
fine tuning Ska, so that the two minima have equal depth for k = 0. The ratio of 
susceptibilities can be computed through the evaluation of the second derivatives of 
the potential at the two minima. Our final estimateis ^X+fX- = 1-7(7). It is in good 
agreement with the predictions of the e-expansionHu (x+/X-rjn 2-0, 2.9, 2.3 for the 
first three orders). However, it disagrees with the lattice resulttffl (x+/X- — 4-1(5)). 

In summary, the approach we presented can provide a quantitative description 
of the universal behaviour near weakly first-order phase transitions. It is based on 
the calculation of a coarse-grained free energy through an exact flow equation. Fixed 
points in the evolution, the appearance of new minima in the potential, and the 
universal properties of the resulting radiatively- induced first-order phase transitions 
can be studied in detail. 



2 Spontaneous nucleation and coarse graining 

2. 1 The necessity of coarse graining 

Near a first-order phase transition the potential (or free energy density) of the sys- 
tem has two separate local minima. As the temperature drops below a critical 
value, the minimum corresponding to the true vacuum (for late times) becomes 
deeper than the one corresponding to the false vacuum. However, the system may 
not adapt immediately to the new equilibrium situation, and we encounter the fa- 
miliar phenomena of supercooling or hysteresis. For sufficiently low temperature 
the phase transition takes place, driven by spontaneous nucleation. For example, 
as vapor is cooled below the critical temperature, local droplets spontaneously form 
and grow until the transition is completed. The inverse evolution proceeds by the 
formation of vapor bubbles in a liquid. The transition in ferromagnets is character- 
ized by the formation of local Weiss domains. The formation of bubbles of the new 
vacuum is similar to a tunnelling process, and typically exponentially suppressed at 
the early stages of the transition. The reason is the barrier between the local min- 
ima. The transition requires the formation of the configuration with lowest action 
on the barrier. The rate includes a Boltzmann factor involving the action of this 
critical configuration. 
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Our present understanding of the phenomenon of nucleation is based largely 
on the work aLLangertJ. His..approach has been appliecLtp relativisticJfield theory 
by Coleman E£l and Callan E£l and extended by Affleck Ell and Linde La to finite- 
temperature quantum field theory. The basic quantity in this approach is the nu- 
cleation rate /, which gives the probability per unit time and volume to nucleate a 
certain region of the stable phase (the true vacuum) within the metastable phase 
(the false vacuum). The calculation of I relies on a semiclassical approximation 
around a dominant saddle-point that is identified with the critical bubble. This 
is a static configuration (usually assumed to be spherically symmetric) within the 
metastable phase whose interior consists of the stable phase. It has a certain radius 
that can be determined from the parameters of the underlying theory. Bubbles 
slightly larger than the critical one expand rapidly, thus converting the metastable 
phase into the stable one. 

The nucleation rate is exponentially suppressed by the action of the critical bub- 
ble. Possible deformations of the critical bubble generate a static pre-exponential 
factor. The leading contribution to it has the form of a ratio of fluctuation determi- 
nants and corresponds to the first-order correction to the semiclassical result. For a 
four-dimensional theory of a real scalar field at temperature T, the nucleation rate 
is given by 



Here V is the effective action of the system for a given configuration of the field <ft 
that acts as the order parameter of the problem. The action of the critical bubble 
is Tb = r [</>b(r)] — r[0], where </>&(r) is the spherically-symmetric bubble configu- 
ration and <j) — corresponds to the false vacuum. The fluctuation determinants 
are evaluated either at (f> = or around 4> = <j>b{r). The prime in the fluctuation 
determinant around the bubble denotes that the three zero eigenvalues of the op- 
erator [5 2 r/<5</> 2 ]0 = 6 have been removed. Their contribution generates the factor 
(rb/27r) 3 / 2 and the volume factor that is absorbed in the definition of / (nucleation 
rate per unit volume). The quantity Eq is the square root of the absolute value of 
the unique negative eigenvalue. 

In field theory, the rescaled free energy density of a system for homogeneous 
configurations is identified with the temperature-dependent effective potential. This 
is often evaluated through some perturbative scheme, such as the loop expansionliil 
In this way, the profile and the action of the bubble are determined through the 
potential. This approach, however, faces three fundamental difficulties: 

a) The effective potential is a convex function of the field and seems inappropriate 
for the the study of tunnelling. 

b) The fluctuation determinants in the expression for the nucleation rate have a 
form completely analogous to the one-loop correction to the potential. The ques- 
tion of double-counting the effect of fluctuations (in the potential and the prefactor) 
must be properly addressed. 

c) The fluctuation determinants in the prefactor are ultraviolet divergent. An ap- 
propriate regularization scheme, consistent with the one used in the calculation of 
the potential, must be employed. 





6 



All the above issues can be resolved through the implemcntion of the notion 
of coarse graining in the formalism E3, in agreement with Langer's philosophy. The 
problem of computing the difference of the effective action between the critical 
bubble and the false vacuum may be divided into three steps: In the first step, 
one only includes fluctuations with momenta larger than a scale k of the order 
of the typical gradients of 4>b(f). For this step one can consider approximately 
constant fields <j> and use a derivative expansion for the resulting coarse-grained free 
energy Tfc [<fi] . The second step searches for the configuration <\>\, (r) which is a saddle 
point of Tk- The quantity Ti, in eq. @ is identified with I^. \<j>b[ — [0] . Finally, the 
remaining fluctuations with momenta smaller than k are evaluated in a saddle-point 
approximation around <f>b (r) . This yields the ratio of fluctuation determinants with 
an ultraviolet cutoff k. Langer's approach corresponds to a one-loop approximation 
around the dominant saddle point for fluctuations with momenta smaller than a 
coarse- graining scale k. 



2.2 Calculation of the nucleation rate 

In the following we review studies of nucleation based on the formalism of the 
effective average action T^, which can be identified with the free energy, rescaled 
by the temperature, at a given coarse-graining scale k. In the simplest case, we 
consider a statistical system with one space-dependent degree of freedom described 
by a real scalar field 4>(x). For example, <j)[x) may correspond to the density for 
the gas/liquid transition, or to a difference in concentrations for chemical phase 
transitions, or to magnetization for the ferromagnetic transition. Our discussion 
also applies to a quantum field theory in thermal quasi-equilibrium. An effective 
three-dimensional description applies for a thermal quantum field theory at scales 
k below the temperature T. We assume that Ffc has been computed (for example 
perturbatively) for ko = T and concentrate here on the three-dimensional (effective) 
theory. 

We compute F^ by solving the flow equation between k and k. For this pur- 
pose we approximate F^ by a standard kinetic term and a general potential U k . 
This corresponds to the first level of the derivative expansion ErlLil. The long-range 
collective fluctuations are not yet important at a short-distance scale k^ 1 = T^ 1 . 
For this reason, we assume here a polynomial potential 

UkoM = \mlj 2 + \ lko f + \x k0 ^ . (7) 

The parameters , 7fe and Xk depend on T . 

We compute the form of the potential Uk at scales k < ko by integrating an 
evolution equation obtained from the exact flow equation for the effective average 
action. Here we employ a mass-like infrared cutoff k for the fluctuations that are, 
incorporated in Tk- The evolution equation for the potential can be written ascSlEll 



v /fc 2 + ^(0)- v /F + ^(O) 
In this entire section primes denote derivatives with respect to 4>. 



(8) 
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The form of Uk changes as the effect of fluctuations with momenta above the 
decreasing scale k is incorporated in the effective couplings of the theory. In general, 
Uk is not convex for non-zero k. It approaches the convex effective potential only in 
the limit k — > 0. In the region relevant for a first-order phase transition, Uk has two 
distinct local minima. The nucleation rate should be computed for k larger than or 
around the scale kf at which Uk starts receiving important contributions from field 
configurations that interpolate between the two minima. This happens when the 
negative curvature at the top of the barrier becomes approximately equal to —k 2 . 
Another consistency check for the above choice of k is the typical length scale of a 
thick-wall critical bubble, which is > 1/k for k > kf. The use of Uk at a non-zero 
value of k resolves the first fundamental difficulty related to the convexity of the 
potential. 

The other two difficulties are overcome as well. In our approach the pre- 
exponential factor in eq. (^) is well-defined and finite, as an ultraviolet cutoff of 
order k must implemented in the calculation of the fluctuation determinants. The 
cutoff must guarantee that fluctuations with characteristic momenta q 2 ?Z k 2 do 
not contribute to the determinants. This is natural, as all fluctuations with typical 
momenta above k are already incorporated in the form of Uk- 

In order to implement the appropriate ultraviolet cutoff ~ k in the fluctuation 
determinant, let us look at the first step of an iterative solution of eq. rtl) 



t r i 1) W-^ 1) (o)= u k M-u k M 



dct[-d 2 + k 2 + U'k'(cj))} det[-d 2 + k 2 + U> k '(0)} 
dct[~d 2 + k 2 + U'l (</>)} det[-d 2 + k 2 + u{{0)] 



.(9) 



For k — > 0, this solution is a regularized one-loop approximation to the effective 
potential. Due to the ratio of determinants, only momentum modes with k 2 < q 2 < 
kg are effectively included in the momentum integrals. The form of the infrared 
cutoff in eq. (||) suggests that we should implement the ultraviolet cutoff for the 
fluctuation determinant in the nucleation rate (^) as 



I = A k exp(-Sfc) 

2 rtet' \-f) 2 + U','(rhJrY)l det \-B 2 + k 2 4- U!'(())l 

(10) 



2tt \2n J 



det [~d 2 + k 2 + E7£'(0 6 (r))] det [~d 2 + U k '(0)] 



We switched the notation to Sk — Tki't'b} — Fk[0] instead of Tb in order to make the 
fc-dependence in the exponential suppression factor explicitly visible. 

As a test of the validity of the approach, the result for the rate I must be 
independent of the coarse-graining scale k, because the latter should be considered 
only as a technical device. In the following we show that this is indeed the case 
when the expansion around the saddle point is convergent and the calculation of the 
nucleation rate reliable. Moreover, the residual k dependence of the rate can be used 
as a measure of the contribution of the next order in the saddle-point expansion. 

The critical bubble configuration (f>b{r) is an SO(3)-invariant solution of the 
classical equations of motion that interpolates between the local maxima of the po- 
tential —Uk{4>)- It can be computed easily with numerical methods. The evaluation 
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Figure 3: Dependence of effective potential, critical bubble and nucleation rate on the coarse 
graining scale k. 

of the fluctuation determinants is more complicated. However, the regular- 
ized expression in (pJL.can be computed through a combination of numerical and 
analytical techniqucsEJEi 

2.3 Two examples 

Sample computations are presented in fig. 3. The first row corresponds to a 
model with an initial potential given by eq. (0) with m 2 . = —0.0433 k 2 , jk = 

—0.0634 kf/ 2 , Afc = 0.1 kg. We first show the evolution of the potential Uk{4>) as 
the scale k is lowered. The solid line corresponds to k/ko = 0.513 while the line with 
longest dashes (that has the smallest barrier height) corresponds to fc//fco = 0.223. 
At the scale kf the negative curvature at the top of the barrier is slightly larger than 
—k 2 and we stop the evolution. The potential and the field have been normalized 
with respect to kf. As A: is lowered from ko to kf, the absolute minimum of the 
potential settles at a non-zero value of <p, while a significant barrier separates it 
from the metastable minimum at 4> = 0. The profile of the critical bubble </>&(r) 
is plotted in the second figure in units of kf for the same sequence of scales. For 
k ~ kf the characteristic length scale of the bubble profile and 1/k are comparable. 
This is expected, because the form of the profile is determined by the barrier of 
the potential, whose curvature is ~ — k 2 at this point. This is an indication that 
we should not proceed to coarse-graining scales below k f . We observe a significant 
variation of the value of the field cf> in the interior of the bubble for different k. 
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Our results for the nucleation rate are presented in the third figure. The hor- 
izontal axis corresponds to kj yJU'^(4>t), i.e. the ratio of the scale k to the square 
root of the positive curvature (equal to the mass of the field) at the absolute min- 
imum of the potential located at <pf Typically, when k crosses below this mass, 
the massive fluctuations of the field start decoupling. The evolution of the convex 
parts of the potential slows down and eventually stops. The dark diamonds give the 
values of the action Sk of the critical bubble. We observe a strong k dependence 
of this quantity. The stars indicate the values of ln(Ak/kj). Again a substan- 
tial decrease with decreasing k is observed. The dark squares give our results for 
— \n(I/kj) = Sk — \n(Ak/kj). It is remarkable that the k dependence of this quan- 
tity almost disappears for k / \/U^((f>t) 1 . The small residual dependence on k 
can be used to estimate the contribution of the next order in the expansion around 
the saddle point. It is reassuring that this contribution is expected to be smaller 
than ln(A k /kj). 

This behaviour confirms our expectation that the nucleation rate should be 
independent of the scale k that we introduced as a calculational tool. It also 
demonstrates that all the configurations plotted in the second figure give equiv- 
alent descriptions of the system, at least for the lower values of k. This indicates 
that the critical bubble should not be associated only with the saddle point of the 
semiclassical approximation, whose action is scale dependent. It is the combina- 
tion of the saddle point and its possible deformations in the thermal bath that has 
physical meaning. 

For smaller values of |m| | the dependence of the nucleation rate on k becomes 
more pronounced. We demonstrate this in the second row in fig. 3 for which 
A feo /(-m^ o ) 1/2 = 0.88 (instead of 0.48 for the first row). The value of A feo is the 

same as before, whereas 7fc = —1.61 • 10~ 3 k^ 2 and kf/k = 0.0421. Higher-loop 
contributions to Ak become important and the expansion around the saddle point 
does not converge any more. There are two clear indications of the breakdown of 
the expansion: 

a) The values of the leading and subleading contributions to the nucleation rate, 
Sk and \n(Ak/kf) respectively, become comparable. 

b) The k dependence of \n(I/kj) is strong and must be canceled by the higher-order 
contributions. 

The discontinuity in the order parameter at the phase transition is approximately 
5 times smaller in the second example than in the first one. As a result, the second 
phase transition can be characterized as weaker. Typically, the breakdown of the 
saddle-point approximation is associated with weak first-order phase transitions. 

In the last figure of each row we also display the values of ln(Ak/kj) (dark 
triangles) predicted by the approximate expression 



Ak ^ irk 
kj 2 



- / r 3 [^(^(r))-^'(0)]dr 
Jo 



1/2 

(11) 



It gives a good approximation to the exact numerical results, especially near kf, 
and can be used for quick checks of the validity of the expansion around the saddle 
point. 
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2-4 Region of validity of homogeneous nucleation theory 

It is useful to obtain some intuition on the behaviour of the nucleation rate by using 
the approximate expression (pi]). We assume that the potential has a form similar 
to eq. (Q) even near kf, i.e. 

UkM) « \ m t f <f> 2 + \lk f ^ + \\k^\ (12) 

(Without loss of generality we take m| > 0.) For systems not very close to the 
endpoint of the first-order critical line, our assumption is supported by the numerical 
data. The scale kf is determined by the relation 



k 2 



US f W =^-m%. (13) 



Through the rescalings r = f/m^, = 2<f)m 2 . ,/jkf> the potential can be 

written as U($) = <p 2 /2 - <p 3 /3 + h4> 4 /!8, with h = 9A fc/ mL/7L. For h » 1 the 
two minima of the potential have approximately equal depth. The action of the 
saddle point can be expressed as 

S kf = --^hS(h), (14) 
9 X kf 

where S(h) must be determined numerically through U (4>). Similarly, the pre- 
exponential factor can be estimated through eq. (0) as 



Ah, 7T / 3 



A 2 {h) = j [&" (* t (f)) - l] f* if, (15) 



with A(h) computed numerically. Finally, the relative importance of the fluctuation 
determinant is given by 

_ in fcyg) ftr 1 jffij *tL ■ rW itL. (16 ) 

Sibj, 8 h V 2/i m fe/ m fc/ 

The ratio i? can be used as an indicator of the validity of the saddle point expansion. 
The latter is valid only for R ^ 1. 

Numerical evaluation of T(h) shows that it diverges for h — > 0, while it becomes 
small for h — > 1. This suggests two cases in which the expansion around the saddle 
point is expected to break down: 

a) For fixed Xk f /mk f , the ratio R becomes larger than 1 for h — > 0. In this limit 
the barrier becomes negligible and the system is close to the spinodal line. 

b) For fixed h, R can be large for sufficiently large \k f /mk f - This is possible 
even for h close to 1, so that the system is far from the spinodal line. This case 



11 



corresponds to weak first-order phase transitions, as can be verified by observing 
that the saddle-point action (ph, the location of the true vacuum 



ffiikj 3 y \ kf 



y/rtSh, 

and the difference in free-energy density between the minima 



(17) 



AC/ =W(/^ (18) 



go to zero in the limit nik f /\k f ~ > for fixed h. This is in agreement with the 
discussion of fig. 3 in the previous subsection. 

The breakdown of homogeneous nucleation theory in both the abpjfe cases is con- 
firmed through the numerical computation of the nucleation rates Eil. 



2.5 Other results 

The approach has been applied to more complicated systems, such as theories of two 
scalar fields. The most interesting feature of the two-scalar models is the presence 
of radiatively induced first-order phase transitions. Such transitions usually take 
place when the expectation value of a certain field generates the mass of another 
through the Higgs mechanism. The fluctuations of the second field can induce the 
appearance of new minima in the potential of the first, resulting in first-order phase 
transitions. We discussed such an example in the first section. The problem of 
double-counting the effect of fluctuations is particularly acute in such situations. 
The introduction of a coarse-graining scale k resolves this problem, by separating 
the high-frequency fluctuations of the system which may be responsible for the 
presence of the second minimum through the Coleman- Weinberg mechanism Ea, 
from the low-frequency ones which are relevant for tunnelling. 

Unfortunately, the expansion around the saddle point does not converge for 
radiatively induced first-order phase transitions. The prefactor resulting from the 
fluctuation determinant of the second field is always comparable to the leading ex- 
ponential term E3. As a result, the saddle-point approximation breaks down and 
the predicted nucleation rate I/kj is strongly k dependent. The above results are 
not surprising. The radiative corrections to the potential and the pre-exponential 
factor have a very similar form of fluctuation determinants. When the radiative 
corrections are large enough to modify the bare potential and generate a new min- 
imum, the pre-exponential factor should be expected to be important also. These 
findings haw important implications for cosmological phase transitions, such as the 
electroweakE3iHJ. 

Finally, the reliability of our approach has been confirmed through compar- 
isons with results obtained through lattice simulations cj, or alternative analytical 
methods in their region of applicability!^. 
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